Does the 1/f frequency-scaling of brain signals 
reflect self-organized critical states ? 
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Many complex systems display self-organized critical states characterized by 1/f frequency scaling of power 
spectra. Global variables such as the electroencephalogram, scale as 1/f, which could be the sign of self- 
organized critical states in neuronal activity. By analyzing simultaneous recordings of global and neuronal 
activities, we confirm the 1 // scaling of global variables for selected frequency bands, but show that neuronal 
activity is not consistent with critical states. We propose a model of 1/f scaling which does not rely on critical 
states, and which is testable experimentally. 
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Self-organized critical states are found for many 
complex systems in nature, from earthquakes to 
avalanches QIQ]. Such systems are characterized by 
scale invariance, which is usually identified as a power- 
law distribution of variables such as event duration or 
the waiting time between events. 1 // noise is usually 
considered as a footprint of such systems |Q] . 1 // fre- 
quency scaling is interesting, because it betrays long- 
lasting correlations in the system, similar to the behav- 
ior near critical points. 

Several lines of evidence point to the existence of 
such critical states in brain activity. Global variables, 
such as the electroencephalogram (EEG) and magne- 
toencephalogram, display frequency scaling close to 
l/f 01- EEG analysis j5| and avalanche analysis of 
local field potentials (LFPs) recorded in vitro |g] pro- 
vided clear evidence for self-organized critical states 
with power-law distributions. There is also evidence for 
critical states from the power-law scaling of interspike 
interval (ISI) distributions computed from retinal, visual 
thalamus and primary visual cortex neurons 0] . In ad- 
dition, model networks of neurons indicate that critical 
states may be associated with frequency scaling consis- 
tent with experiments 0. However, these are indepen- 
dent evidences from different preparations and the link 
between 1 // frequency scaling of global variables and 
the existence of critical states in neural activity has not 
been firmly established. Moreover, 1 // spectra are not 
necessarily associated with critical states |9], so it is not 
clear if the intact and functioning brain operates in a 



way similar to critical states. 

To attempt answering these questions, we first inves- 
tigated if 1 // frequency scaling is present in global vari- 
ables recorded close to the underlying neuronal current 
sources in vivo. We analyzed cortical LFPs which were 
recorded within cerebral cortex using bipolar extracellu- 
lar high-impedance microelectrodes 11(11 . Bipolar LFP 
recordings sample relatively localized populations of 
neurons, as these signals can be very different for elec- 
trodes separated by 1 mm apart 1 10]. This stands in con- 
trast with the EEG, which samples much larger popula- 
tions of neurons II ill and is recorded from the surface 
of the scalp using millimeter-scale electrodes. LFPs are 
subject to much less filtering compared to EEG, be- 
cause EEG signals must diffuse through various me- 
dia, such as cerebrospinal fluid, dura matter, cranium, 
muscle and skin. Thus, finding l/f frequency scal- 
ing of bipolar LFPs would be a much stronger evidence 
that this scaling reflects neuronal activities, as these sig- 
nals are directly recorded from within the neuronal tis- 
sue. Moreover, in order to distinguish state-dependent 
scaling properties, we have compared recordings during 
wakefulness and slow-wave sleep in the same experi- 
ments. 

Bipolar LFPs from cat parietal association cortex 
show the classic landmarks of EEG signals in these 
states Il2ll . namely during waking, LFPs are of low 
amplitude and very irregular (Fig. ^ top trace), and 
are dominated by beta frequencies (around 20 Hz). 
This pattern is also called "desynchronized" activity, 
and is typically seen during aroused states in the hu- 
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FIG. 1: Frequency scaling of local field potentials from cat 
parietal cortex. Top traces: LFPs recorded in cat parietal cor- 
tex during wake and slow-wave sleep (SWS) states. Bottom: 
Power spectral density of LFPs, calculated from 55 sec sam- 
pled at 300 Hz (150 Hz 4th-order low-pass filter), and rep- 
resented in log-log scale (dashed lines represent l// a scal- 
ing). During waking (black), the frequency band below 20 Hz 
scales approximately as 1// (*: peak at 20 Hz beta frequency), 
whereas the frequency band between 20 and 65 Hz scales 
approximately as l// 3 . During slow-wave sleep (gray; dis- 
placed upwards), the power in the slow frequency band is in- 
creased, and the 1// scaling is no longer visible, but the 1 // 3 
scaling at high frequencies remains unaffected. PSDs were 
calculated over successive epochs of 32 sec, which were aver- 
aged over a total period of 200 sec for Wake and 500 sec for 
SWS. 

man EEG fiUl . During slow- wave sleep, LFPs dis- 
play high-amplitude slow-wave activity (Fig. ^ mid- 
dle trace), similar to the "delta waves" of human sleep 
EEG |11]. The power spectral density (PSD) calculated 
from these LFPs typically shows a broad-band struc- 
ture. During wakefulness, the PSD shows two different 
scaling regions, according to the frequency band. For 
low frequencies (between 1 and 20 Hz), the PSD scales 
approximately as 1 //, whereas for higher frequencies 



(between 20 and 65 Hz), the scaling is approximately 
of 1 If (Fig. G] black PSD). During slow -wave sleep, 
the additional power at slow frequencies masks the 1 // 
scaling, but the same 1 // 3 scaling is present in the high- 
frequency band (Fig.^ g ra y PSD). The same behavior 
was observed for other electrodes in the same experi- 
ment, and in three other animals (not shown). Thus, 
these results confirm that the 1 // frequency scaling re- 
ported in the EEG Ql is also present in bipolar LFPs 
from cat association cortex, but only during waking and 
for specific frequency bands. 

To investigate whether this l/f scaling is associated 
with self-organized critical states, we first analyzed the 
ISI distributions from neurons recorded in cat parietal 
cortex. Unit activity was recorded simultaneously with 
LFPs at 8 locations separated by 1 mm The dis- 
tribution of ISIs was computed for individual neurons, 
and were represented in log-linear scale (Fig. |2j log- 
log scale in insets). For both wakefulness and slow- 
wave sleep (Fig.|2j\ and B), the distributions showed no 
evidence for power law behavior. During waking, the 
ISI distributions were close to exponentially-distributed 
ISIs, as generated by Poisson stochastic process with 
same statistics as the neurons analyzed (Fig. [2] Poisson). 
For 22 neurons recorded during the wake state, the Pear- 
son coefficient was of 0.91 ±0.13 for exponential dis- 
tribution fits, and of 0.86 ± 0. 16 for power-law distribu- 
tion fits. Taking only the subset of 7 neurons with more 
than 2000 spikes, the fit was nearly perfect for exponen- 
tial distributions (Pearson coefficient of 0.999 ±0.001). 
However, during slow-wave sleep, there was a marked 
difference between the experimental ISI and the corre- 
sponding Poisson process (Fig. |2|3). In this state, neu- 
rons tended to produce long periods of silences, which 
are related to EEG slow waves and which is 

visible as a prominent tail of the distribution for large 
ISIs. This tail was well fit by a Poisson process of low 
rate (Fig. |2j3, dashed line). 

To further check for criticality, we have performed 
an avalanche analysis by taking into account the col- 
lective information from the multisite recordings. We 
used the same method as for ref. |6], which amounts 
to detect clusters of contiguous events separated by si- 
lences, by binning the system in time windows of 1 ms 
to 16 ms |6J]. As there was no evidence for any recogniz- 
able event in LFPs which could be taken as avalanche 
(see Fig. we used the spike times among the en- 
semble of simultaneously recorded neurons. The dis- 
tribution of avalanche size does not follow power-law 
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scaling (Fig.|2jZ!, black), but is closer to an exponential 
distribution as predicted by Poisson processes (Fig. |2jT, 
gray). This analysis therefore confirms the absence of 
avalanche dynamics in this system filill . 

To explain the 1 // scaling of LFPs, we attempted to 
reconstruct LFPs from unit activity. Unit activity is dis- 
played in Fig.|5](top) for the same experiment as that of 
Fig-H Because LFPs are generated primarily by synap- 
tic currents in neurons llll Il3ll . and because synaptic 
currents are very well modeled by simple exponential 
relaxation processes fl4Tl . we modeled the synaptic cur- 
rent from the following convolution Hill : 

C(t) = J~ D(t') exp[—(t — t')/x s )] dt' , (1) 

where C(t) is the synaptic current and D(t) is the "drive" 
signal which consisted in the experimentally-recorded 
spike trains. The PSD of the synaptic current is then 
given by: 

S((0) = \C{(0)\ = 1+& 2 Z 2 ■ ( 2 ) 

The PSD of synaptic currents reconstructed from 
experimentally-recorded spikes showed an approximate 
Lorentzian behavior (l/f 2 scaling) during wakefulness 
(Fig. [3] Wake), as expected from the exponential na- 
ture of synaptic events. During slow-wave sleep, there 
was more power for slow frequencies, but the l/f 2 scal- 
ing at high frequencies was still present (Fig. [3] SWS). 
The Lorentzian form of the PSD in Fig.[3](Wake) shows 
that in the waking state, |D(co)| 2 is approximately con- 
stant, therefore the drive D(t) is statistically equivalent 
to a white noise process, consistent with the apparent 
Poisson statistics of spikes identified in Fig.|3](see also 
refs. fl7il for similar findings in awake monkeys). Dur- 
ing slow-wave sleep, however, the deviation from the 
Lorentzian suggests that D(t ) is a stochastic process sta- 
tistically different from white noise, and contains in ad- 
dition increased power at low frequencies, also consis- 
tent with the analysis of Fig- E] 

This model, however, does not yield PSD consistent 
with the l/f and l// 3 scaling of LFPs shown in Fig.^ 
Interestingly, the scaling of this model is in 1//° or 
1 // 2 for the same frequency bands that displayed l/f 
or l// 3 in LFPs, respectively. Using a similar convolu- 
tion equation to model the PSD of LFPs 

LFP{t) = ( C{t')F{t-t')dt' , (3) 



where C(t) is the synaptic current source and F(t) is 
a function representing a filter. As above, the PSD is 
given by 

P((o) = \LFP(a)\ 2 = |C(co)| 2 \F(o))\ 2 . (4) 

In this model, the frequency scaling of the PSD of 
both wakefulness and slow-wave sleep LFPs in Fig. [2 
can be explained by assuming that the filter scales as 
l/f, or equivalently that ^(d))! 2 ~ 1/(0. In other 
words, this model can explain qualitatively the 1 // and 
l// 3 scaling of LFPs under the condition that neuronal 
current sources are subject to an 1 // filter. Such a fil- 
ter is most likely due to the filtering of extracellular 
currents through the tissue, before it reaches the elec- 
trode Hii . 

Finally, we provide an intuitive justification for this 
predicted 1 // filter, as well as possible ways to test it 
experimentally. The l/f filtering of extracellular me- 
dia can be justified intuitively by considering the com- 
plex structure of such media, and in particular its spa- 
tial irregularity. Extracellular space consists of a com- 
plex arrangement of cellular processes of various size 
and irregular shape, while the extracellular fluid repre- 
sents only a few percent of the available space 1 1 911 . The 
effect of a current source in such media will be a com- 
bination of resistive effects, due to the flow of current 
in the conductive fluids, and capacitive effects, due to 
the high density of membranes (for a theoretical treat- 
ment see refs. JTill '). Such a complex arrangement of 
resistors and capacitors with random values is known 
to produce an 1 // filter, as found for inhomogeneous 
materials lEoll . Although such materials are different 
from the structure of biological media, it is plausible 
that similar considerations may explain the 1 // filtering 
predicted here. Linear arrangements of RC circuits with 
random values (RC line) also generates l/f noise I2II1 . 
Superposition of a large number of exponential relax- 
ation processes with different relaxation rates can also 
generate 1 //scaling I22H23I1 . Understanding of the l/f 
filtering by extracellular media based on plausible bio- 
physical models is presently under investigation. The 
predicted l/f filter could also be tested experimentally 
by injecting white noise currents (of amplitude compa- 
rable to neuronal current sources) in extracellular space, 
and measuring the resulting field potential at some dis- 
tance from the injection site. This measured LFP should 
scale as 1 //. 

In conclusion, we have shown that the PSD of bipo- 
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lar LFPs from cat parietal cortex displays several scal- 
ing regions, as 1// or l// 3 depending on the frequency 
band and behavioral state. By analyzing neuronal unit 
activity from the same experiments, we did not see ev- 
idence that this l/f scaling is associated with critical 
states. Neither ISI distributions nor avalanche size dis- 
tributions display power-law scaling, but are rather con- 
sistent with Poisson processes. We provided an alterna- 
tive explanation for l/f frequency scaling which does 
not rely on critical states, but rather stems from the fil- 
tering properties of extracellular media. We gave an 
intuitive explanation for a possible physical origin of 
such l/f filtering, as well as a way to test it experi- 
mentally. These results may appear to contradict with 
previous evidence for critical states in vitro J(J or in the 
early visual system in vivo 0|. However, the absence of 
critical states reported here may instead reflect funda- 
mental differences between association cortex and other 
structures more directly related to sensory inputs. Fu- 
ture work should clarify why different structures show 
different scaling, and what implications it may have for 
brain dynamics and coding. 
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FIG. 2: Absence of power law distributions in neuronal ac- 
tivity. The logarithm of the distribution of interspike inter- 
vals (ISI) during waking (Wake, A, 1951 spikes) and slow- 
wave sleep (SWS, B, 15997 spikes) is plotted as a function 
of ISI length, or log ISI length (insets). A poisson process 
of the same rate and statistics is displayed in A (Poisson; 
gray curve displaced upwards for clarity). The exponential 
ISI distribution predicted by Poisson processes of equivalent 
rates is shown as straight lines (smooth curve in inset). The 
dotted line in B indicates a Poisson process with lower rate 
which fits the tail of the ISI distribution in SWS. C. Avalanche 
analysis realized by taking into account the statistics from all 
simultaneously-recorded cells in Wake. The distribution of 
avalanche sizes scales exponentially (black curves), similar to 
the same analysis performed on a Poisson process with same 
statistics (gray curves). 
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FIG. 3: Frequency scaling of synaptic currents reconstructed 
from spike times. Top traces: raster plot of spiking times 
of 8 multi-unit recordings in cat cortex during wakefulness 
(same experiment as in Fig. Q data from ref. IUm ). Mid- 
dle trace: total synaptic current obtained by convolving the 
spike times with exponential relaxation processes (t s =10 ms). 
Bottom: PSD of synaptic currents for wake (black) and slow- 
wave sleep (SWS; PSD in gray displaced upwards for clarity); 
dashed lines represent 1 / f a scaling. 



